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Abstract 

We consider the one-dimensional delta-interacting electron gas in 
the case of infinite repulsion. We use determinant representations 
to study the long time, large distance asymptotics of correlation 
functions of local fields in the gas phase. We derive differential 
equations which drive the correlation functions. Using a related 
Riemann-Hilbert problem we obtain formulae for the asymptotics 
of the correlation functions, which are valid at all finite tem- 
peratures. At low temperatures these formulae lead to explicit 
asymptotic expressions for the correlation functions, which de- 
scribe power law behavior and exponential decay as functions of 
temperature, magnetic field and chemical potential. 
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1 Introduction 

We consider the electron gas with delta interaction in one space and one 
time dimension. It can be described by canonical Fermi fields ip a {x) with 
canonical equal time anticommutation relations 

M*H(y) + 4(y)M*) = tfs{x - y ) . (i) 

The spin index a runs through two values, a =]■,[■ The Hamiltonian of the 
model is 

H = dn{d x $ld x $ a + c:(fj,^ . (2) 

Here [i is the chemical potential, and B is the magnetic field. a z is a Pauli 
matrix, and c is the coupling constant. In this paper we shall consider the 
case of infinite repulsion c = oo. 

The model can be solved exactly Its thermodynamics was described 
in 0. The pressure 

T roo 

P = — dk\n(l + e^ +B - k ^ T + e ^- B - k ^ T ) (3) 

serves as thermodynamic potential. T is the temperature. Note that the ex- 
pression for P, eq. @, is formally the same as for free fermions with effective 
chemical potential fi + T ln(2 cosh(5/T)). Thus, there are two different zero 
temperature phases depending on whether limT^o+(A i + ^l n (2 cosh(i?/T)) = 
/i + \B\ is positive or negative. For fi + \B\ > the density D = dP/dfi 
has a positive limit as T goes to zero. For \i + \B\ < the density at zero 
temperature vanishes. This is the phase we are interested in. We call it the 
gas phase. We shall assume in the following that /i + Tln(2 cosh(i?/T)) < 0. 

The number of up-spin particles and the number of down-spin particles 
are separately conserved. The densities of up-spin electrons and of 
down-spin electrons are obtained as 

OP J_ [°° g— B/T 

~ d(fi-B) 2n J-oo 2cosh(B/T) +e( k2 ~^/ T ' [ ) 

dP 1 r°° e B/T 

° l = d(fi + B) = 27 loo dk 2cosh( J B/T) + e( fc2 -^ ' (5) 
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It is obvious from these expressions that they satisfy the relation 

D l (B) = D^-B) . (6) 

The density D of the gas is 

Note that the density vanishes exponentially at small temperatures, 



D 



^e^ T B = , 

if (8) 

1 T e ^+\B\)/T B , Q 

V Z V 7T ' 



At very low temperatures density and pressure are related by the ideal gas 
law, 

P = DT . (9) 

For larger temperatures there are corrections to the ideal gas law. 

In this article we study the two-point correlation functions G^(x, t), a =|, 
I, of local fields. We shall follow the approach of the book [|3| to correlation 
functions of solvable models. The definitions of G^(x,t) are given in the 
next section, where we recall a determinant representation for G^(x, t), ||, |], 
which is the starting point of our considerations. The correlation functions 
G^(x,t) are not elementary functions. They are determined by a pair of 
coupled nonlinear partial differential equations || which is derived in section 
3. In section 4 we shall formulate a corresponding Riemann-Hilbert problem 
M- H Hi • ft wn l fi x t ne solution of the differential equations, which is relevant 
for the calculation of the correlation functions. Furthermore, the asymptotic 
solution of the Riemann-Hilbert problem for large x and t (x/t fixed) can 
be obtained by reducing it to a known form [|1(], [□], [12], [13], [14], [TJ]. This 



is done in section 5. Using the asymptotic solution of the Riemann-Hilbert 
problem and the differential equations we shall obtain the asymptotics of 
the correlation functions G^(x,t) in terms of a contour integral in section 6. 
In section 7 we shall present explicit asymptotic expressions for G„(x,t) at 
low temperatures. The asymptotic expressions for the correlation functions 
G^(x,t) at low temperatures are products of correlation functions of free 
fermions and corrections caused by the interaction. The leading factors of 
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the correction terms describe exponential decay. The next to leading factors 
are power law corrections. The leading factor of the asymptotics of the spin- 
up correlation function Gf(x,t), for instance, is 

The factor g^O-- 8 )^ 2 /^ j s the exponential part of the free fermionic cor- 
relation function. The factor e~ xL>l appears due to the interaction. The 
explicit form of the power law correction as a function of temperature, chem- 
ical potential and magnetic field will be given in section 7 (see eqs. ( |133| ), 
( |132|) , ( |135|) , (|128|) ). The occurrence of the density of down-spin electrons in 
equation ([H]) has a clear physical interpretation. Due to the Pauli principle 
and the locality of the interaction, up-spin electrons are scattered only by 
down-spin electrons. D± appears as the reciprocal of the correlation length 



in eq. (10) . The higher the density of down-spin electrons, the more often an 



up-spin electron is scattered and the smaller is the correlation length. 



2 A determinant representation for correla- 
tion functions 

We shall consider the following two-point correlation functions of local fields, 

- tr(e-^ a (M)^t (0 ,o)) 
^ a \ x : 1 ) - tr(e~ H / T ) ' 

tr (e-«/^,t(,, t)fe (0,0)) 
^aK^ L ) - ti(e- H / T ) ' 

where a =|, j. These correlation functions depend not only on space and 
time variables x and t, but also on the temperature T, the chemical potential 
fi (or the density D) and the magnetic field B. Due to the invariance of the 
Hamiltonian (0) under the transformation \ v^|, B # —B, we have the 
identity 

Gf(x,t) = Gf(x,t)\^_ B . (13) 

Hence, we can restrict our attention in the following to the correlation func- 
tions Gt(x, t). These correlation functions were represented by means of 
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determinants of integral operators [||, pfl. Now we shall recall these formulae: 

G <^> = -^L d "T^k"» A < I + '' v ) • (14) 

Cf^') = 6 Z B jyvF(7,v)B-det(i + 7 v) . (15) 
It will take us some time to define notations. First of all 



e 



at 



F(rf,ri) = l+ - + — (16) 

7 — e ir > 7 — e ir ? 

and 

1 = l + e 2B ' T . (17) 

V is an integral operator with kernel V(k,p). It can be represented in the 
"standard form" (see page 316 of the book ||) 

lV ^ P ) = e+{k)e - {P }- eMe ~ {k) . (18) 

k — p 

Here 



e_(*) = JT^e^y* , (19) 



e+ (jfe) = I W^5 e -r(fc)/3 1 (1 _ C0S ^yW E ( k } + sin (^)| ; (20) 

Wj " 2cosh( J B/T)+e( fc2 -^ ' 
r(Jfe) = ikH-ikx (22) 

and 

/•oo e — r(p) 

E(k)=p.v. dp- . (23) 

J-oo 7r(j) — fcj 

The integral operator (|TS|) closely reminds an integral operator which appears 
in the theory of correlation functions of the impenetrable Bose gas (cf. page 
346 of 0). The inverse of (/ + jV) can be defined in the following way, 

(/ + jV) (l-jR)=I . (24) 
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The resolvent R is an integral operator with kernel R(k,p). It can be repre- 
sented in the "standard form" , 

M = fMPlzIMfM (25) 

K p 

where the functions f±{k) can be defined as solutions of the integral equations 

/oo 
dpjV(k,p)f ± (p) = e ± (k) . (26) 
-oo 

This theorem is proven on page 318 of the book ||. We should mention that 
both kernels V(k,p) and R(k,p) are symmetric, 

V(k,p) = V(p,k) , R(k,p) = R(p,k) . (27) 

At this point it is convenient to define potentials B a t and C a bj 

/oo roc 
dke a (k)f b (k) , C ab = dkke a (k)f b (k) . (28) 
-oo J ~ oo 

Here both indices a and b run through two values, a, b = ±. For example, 
the factor B in (pl|) is defined as 

/oo 
dke4k)f.(k) . (29) 
-oo 

In order to define b ++ we need the function 

G(x,t) = -£- dke~ T ^ = -^e lx/it . (30) 



The factor b ++ in ([L4]) is defined as 

/oo 
dke + (k)f + (k)-(l-cos(rj))G(x,t) . 
-oo 

(31) 

In such a way we accomplished the description of the determinant represen- 
tation of the correlation functions (|l4|) , (|l5|). It reminds the determinant 
representation for correlation functions of the impenetrable Bose gas (see 
page 346 of the book ||). Let us also mention that the symmetry of V(k,p) 
leads to B. = B_*. 
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Let us discuss the factor 1— cos(^) in the denominator of QUI). Let us show 
that there is no singularity at r\ = 0. Indeed, e+(k), eq. (pOj), is proportional 
to ?7 as 77 — > 0. Hence 6 ++ ~ r/ 2 , and there is no singularity. Let us change 
the integration variable to z = e 171 . Then 

G+(x,t) = ——j>dzj^-^b ++ (z)det(l + >yV)(z) , (32) 

p -\t(jj,-B) , T?{ z \ „ N 

Gf(z,t) = 4?ri7 f dz -^-B^(z)det (/ + -yV) (z) . (33) 

Here the contour of integration is the unit circle, and 

F(z) = l + ^ + ^— . (34) 
7 — z syz — 1 

3 Differential equations 

Starting from a determinant representation for a quantum correlation func- 
tion it is possible to derive differential equations which drive this correlation 
function (see || and chapter XIV of the book 0). In our specific case calcu- 
lations similar to the ones presented on pages 345 - 353 of the book lead 
us to the following differential equations, 

<%++ d 2 b ++ 



d 2 B 



(35) 

+ 2b ++ B 2 _ _ . 



dt dx 2 
We should notice that a change of variables, 

x = -2x , t = 2t , (36) 

leads to coincidence of our system of eqs. fl3"5]) with the system (5.47) on 



page 344 of the book. The system of equations (^) is called separated 
nonlinear Schrodinger equation. It is completely integrable. It has many 
solutions. The modern way to fix a solution of the system ([35|) is to formulate 
a corresponding Riemann-Hilbert problem (see chapter XV of the book ||). 
We shall formulate this Riemann-Hilbert problem in the next section. We 
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will be able to solve the Riemann-Hilbert problem asymptotically as x — > oo 
and t — > oo (x/t fixed). This will give us asymptotic expressions for b ++ 

and B The solution of the Riemann-Hilbert problem will also give us 

asymptotic expressions for the other potentials B ab and C a b (f28|). 

Let us now explain how to express det( I+jV) in terms of these potentials. 
The logarithmic derivatives of the determinant can be expressed in terms of 
B and C, 

d x In det (/ + 7^) = i£+_ = LB_+ , 

(37) 

d t In det (i + jV) = -i(C+- + C_+ + (1 - cos(r/))Gfi__) . 

We shall also mention several identities, which will be useful for further 
calculations: 

d x B + _ = -ifl__6 ++ , (38) 
= id x B__ - B„_B + _ , (39) 
C++ = -id x B ++ - 2(1 -cos(ri))GB + _+B + _B ++ , (40) 
C+_-C_+ = B 2 + _-B ++ B^ . (41) 

In the next section we shall introduce a Riemann-Hilbert problem and explain 
how one can express the potentials B ab and C ab in terms of its solution x(k)- 



4 Riemann-Hilbert problem 

The modern way to solve completely integrable differential equations is by 
means of an equivalent Riemann-Hilbert problem 0, || ||. Solving a Rie- 
mann-Hilbert problem amounts to constructing a piecewise analytic matrix- 
function. Let us formulate the Riemann-Hilbert problem connected to the 
separated nonlinear Schrodinger equation (PS]). Let us introduce a 2 x 2- 
matrix x(^)- The matrix x{k) should be analytic as a function of k in the 
upper half plane, Im(fc) > 0, and in the lower half plane, lm(k) < 0. In 
general, it has different limits x+{k) as & approaches the real axis from above 
and X-(^) as k approaches the real axis from below. These limits are related 
by means of a conjugation matrix G(k), 

x _(k) = x +(k)G(k) for Im(Jfe) = . (42) 
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The conjugation matrix G(k) is of dimension two. It is defined only on the 
real axis. The matrix should also satisfy the normalization condition 

&xW=I= g J . (43) 

The problem is to find such a matrix m order to relate x{k) to the 

integral operator flTSP we shall use a conjugation matrix of the special 
form 

A 0\ / -e_(jfc)e + (ife) e= (Jfe) \ 

G(k) = +2vri W +W +V 7 . (44) 

\0 V ^ -e 2 _(fc) e + (fc)e_(lfc) J 

The Riemann-Hilbert problem with this form of conjugation matrix was stud- 
ied in detail in chapter XV of the book ||. Let us recall some fundamental 
facts: 

One can prove that 

The functions f±(k) were defined by the integral equations (p6|). On the 
other hand, the Riemann-Hilbert problem can be solved by means of these 
functions 



(I 0\ [°° dp 1 



1/ J-ooj) — k 



e-(p)/+(p) -e+(p)/+(p) 
e-(p)/-(p) -e+(p)/_(p) 



(46) 



In order to express the potentials 5 a t and C a b in terms of x(/c), let us study 
the asymptotics of as — * oo, 



Comparing fllBQ, (KTT) and (BBI) we obtain 



In such a way we have formulated the Riemann-Hilbert problem and related 
it to the integral operator V and to the correlation functions Gt(x,t) and 
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Gj(x,t). Next we shall solve the Riemann-Hilbert problem at large space 
time separation x — > oo, t — > oo, x/t fixed. Then we shall calculate the 
potentials B ab and C a b from pT|). The knowledge of the potentials will per- 
mit us to evaluate det(J + jV), see (j37[). This will give us an asymptotic 
expression for the correlation functions (|14]), (|15|) - We also want to note that 

detG(k) = detx(k) = 1 . (49) 



5 Asymptotic solution of the Riemann-Hil- 
bert problem 

In order to solve the Riemann-Hilbert problem asymptotically ]7|, |K], [11, [12, 
|l~3| , [L4], [L5|] , we shall reformulate it twice. This will bring it into the canonical 



form, which was solved by S. V. Manakov ||T0| , pTJ] . 

First reformulation: We want to simplify the conjugation matrix G(k). 
Let us define a matrix 



Xo{k) 



(x 



— (1 — cos(?7)) 



-oo 2n(p — k) 







(50) 



Here r(p) = ip 2 t — ipx. The matrix Xo(^) is analytic for Im(/c) > and for 
Im(A;) < 0. At large k it approaches the unit matrix as a Taylor series in 
1/k. Now, instead of the unknown matrix function x{k), let us introduce a 
new unknown matrix function x(k), 



x (k) = x{k)xo{k) ■ 
Let us formulate the Riemann-Hilbert problem for x(k). 

(a) x(k) is analytic for Im(fc) > and for Im(fc) < 0. 

(b) As k 



(51) 



oo 



k k 2 



k 3 



(52) 



[c) x(k) has a discontinuity across the real axis, which is described by a 
conjugation matrix G(k), 



X_(k) = x+{k)G{k) for Im(fc) = 



(53) 
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Let us calculate 

I G n (k) G 12 (k) \ 

G(k)=\ 11 U UxKWteW)" 1 • (54) 

y G 2 i(«) G 22 (fc) y 

Substituting (03) and (M) into (pi) we obtain 



On (Jfe) = 1 + r&(k)(e- ir > - 1) , Gi 2 (fc) = -i(l - cos(r]))e- T ^(l - 7 0(A;)) , 
G 2 i(fc) = -2i 7 ^(A;)e r W , G 22 (£;) = 1 + 7 ^(A;)(e i " - 1) , 

(55) 

where 

rt(k) ~ 2cosh(g/T) f 56 ) 
7W "2cosh(S/T)+e( fc2 -^ • (56) 

Now the conjugation matrix is an elementary function of k, det G(k) = 1 and 
detx(k) = 1. 

Using ( |50"1) and (0) we can also evaluate the coefficients of the expansion 
Q52|) of x(^) for large fc, 



and 



X2 

Let us recall that 



C_+ G ++ + (l-cos(?7))(ia K G + G£ + _) 
I -G__ C+_ + (1 - cos(r/))G5__ 



(58) 



1 yOG p — i7r/4 

G(x,t) = — dke- r{ ^ = -^e lx2 l u . (59) 



W=exp|^y~ ^ln(l + 7 ^)(e- i "-l))| . (60) 



In order to prepare for the second step let us introduce a function a(k), 

dp 

-oo p 

This function solves the scalar Riemann-Hilbert problem 

a _(k) = a + {k) (l + ^(k){e"^ - 1)) . (61) 
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The conjugation function here is the element Gn(k) of the conjugation matrix 
G(k), eq. (|55|) . Later we shall need the coefficients a±, a 2 of the expansion 
of ln(a(fc)) for large k, 

Ha{k)) = T + 1 + ° (p) • (62) 

Second reformulation: We want to bring the Riemann-Hilbert problem 
to the canonical form. Instead of let us define another 2 x 2-matrix 

$(k) = x(k)e- aZlna{k) . (63) 

This matrix function is analytic for lm(k) > and lm(k) < 0. It approaches 
the unit matrix at large k, 

m =I + ^ + ^ + 0(±) . (64) 



k k 2 \k 3 

It has a discontinuity on the real axis, which can be described by a conjuga- 
tion matrix G$(k), 

$_(jfc) = $ + (jfe)G$(fc) . (65) 

Here 

The matrix elements of G$(k) are 

(G r $)n = 1 , (G$)i2 = ct + (£;)«_ (£;)Gi 2 (A: ) , 

(G*) 2 i = f, 2 ! (fc L , (G*h2 = G 11 (k)G 22 (k) . 
a-(k)a + {k) 

We also note that 

det G^k) = det$(fc) = 1 . (68) 
Let us introduce the notations 

p(Jfe) = -i(l-cos(»7))(l-7i?(ife))a + (ife)a_(A;) , (69) 
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Then we can write the conjugation matrix Gq>(k) in the canonical form 

/ 



G*(ife) 



1 p{k)e~ T{k) 
q(k)e T ^ l+p{k)q(k) 



(71) 



Here r(k) = ik 2 t — ikx. In the following the stationary point 

/ x 

Kn — — 

2t 

of the phase r(k), df-riko) = 0, will play an important role. 
Let us express the potentials B a b in terms of $i, eq. (|6^). 



(72) 



- ot x b ++ 

-B__ B + _ + ai 



(73) 



Here 



i = — — : f°° dp In (l + 7 ^(p)(e- i " - 1 

/7T1 J-oo v 



(74) 



is the first coefficient of the expansion ( |6^) of ln(a(fc)) for large k. Using ([37]) 
and (|73|) we obtain 

d x In det (/ + jV) = LB+_=ifl_ + = -i(($i)n + on) . (75) 
The time derivative of In det can be expressed in terms of $ 2; ec i- 

©, 

<9 t In det (/ + 7 V") = -i(C+_ + C_+ + (1 - cos(r7))Gfi__) 

= i {tr(or z $ 2 ) + 2« 2 } , (76) 
where a 2 is the second coefficient in the expansion ([62]), 



«2 



2vri 



dppln (l + 7^(p)(e~ 



i'/ 



(77) 



This is the end of the second reformulation. 

The asymptotic solution of the Riemann-Hilbert problem with conjuga- 
tion matrix G$(k), eq. (|7T|), is known. It is given by Manakov's Ansatz (see 
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page 452 of the book ||). In order to describe Manakov's Ansatz we first 
solve a scalar Riemann-Hilbert problem for a function 5(k), 

5 + (k) = 54k) (1 + p(k)q(k)6(k - k)) . (78) 

Here 9 is the Heavyside step function, 



The solution is 



f 1 /■* 
W=exp|- j /_ 



ln(l+p(u) ? (u)U . (80) 



u — k 

Now we can define functions 7 p (fc) and I q {k), 

I p (k) = — -^r 5 + (u)54u)p(u)e-^ , (81) 
2m J-oo u — k 

Mk) = 1 r du q ^ e{U) . (82) 
27ri J-oo u — k 8 + (u)5-(u) 

Using these functions we can write the solution of the Riemann-Hilbert prob- 
lem (ID, ([n]) as 



/ 

^ ta *(*)+0(H) . (83) 



V 



1 -J p (A;) 
-I«(Jfe) 1 



e 



This solution is an asymptotic solution which is valid only as x — » oo, £ — > oo 
for fixed ratio x/t. 

6 Asymptotics of the correlation functions 

We solved the Riemann-Hilbert problem. Now we can calculate the logarith- 
mic derivatives of the determinant, 

d x lndet (1 + 7 V) = LB+_ = LB_ + = -i(($i) n + at) , (84) 
d t lndet (l + 7 v) = i {tr(a"$ 2 ) + 2a 2 } . (85) 
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In order to obtain $1 and $2 we have to decompose our asymptotic solution 
C |83f) into a Taylor series in 1/k fl64|) . Then 

- + . (86) 

For a.\ see (ff4|). In order to define 8\ let us recall that 



MW) = ^r/*° ln (1-2(1 - oos(t7))7^(p) (1-7^))) , (87) 

k = x/2t. Here we used (|S9|) and (|30|). At large k, ln(5(k)) can be decom- 
posed into a Taylor series in 1/k, 

In(tf(*)) = | + | + 0(1) . (88) 



So 

<Ji = -— / dpln(l-2(l-cos(7/))7i?(p)(l-77?(p))) , (89) 
2ix\ J -00 



1 f k o 

6 2 = :/ dppln(l-2(l-cos(r7))7i?(p)(l-7i?(p))) . (90) 



2vri 

Substitution of ai, eq. (|74D, and #1 into (|84]) gives 
5+_ = -— /°° rfpsign(p - jfc ) In fl + 70(p) ( e -w0>-*o) _ ^ ? ( 91 ) 

Z7T1 J-oo ^ ^ ' ' 

where k = x/2t and sign(fc) is the sign function, 

sign(fc) = I if . (92) 

1-1 k < 

Using (|34|) we obtain 

a c lndet(/ + 7V r ) = 

- [°° dp sign(p - fc ) In f 1 + T 0(p) ( e -™(P- fc ") - l) ) . (93) 

27T V-oo v v 7 7 
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We shall calculate the time derivative in a similar way, 

tr(<r*$ 2 ) = 26 2 , (94) 
see (^) and (|9TiD. From flgBj ) we have 

a t lndet(/ + 7 F) = 2i(5 2 + a 2 ) , (95) 
see (^) and ©. Finally, 
5 t lndet(/ + 7l/) = 

- y"°° dppsign(p - fc ) In (l + y&{p) (V™(P- fc °) - l)) . (96) 

Now we have calculated space fl93|) and time derivative of In det(I + r yV). Let 
us integrate, 

lndet(J + jV) = — f°° dp\x- 2pt\ In (l + 7 0(p) ( e ~ i,?sign{p - fco) - l)) . 

(97) 

This is the leading term of the asymptotics. Let us recall that 

im = 2cosh ^) (98) 

7 W 2cosh( J B/T)+e( fc2 ^)/T " W 

In order to evaluate the correlation functions (|I~4|), ( |i~5|) we also need b ++ 
and B 

6 ++ =($i)i2 , B— = -($0 2 i , (99) 

see (|75[) . Using the explicit estimate for the error (see page 453 of the book 
[0]) of the asymptotic solution eq. (83), of the Riemann-Hilbert problem 
we find 

, |B_|~i , (100) 

which shows that 6 ++ and B decay with time. In order to get a more 

precise estimate of the asymptotics, we shall turn to the differential equations 
!5T). We shall use some information from the book H. So we shall rewrite 
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the differential equations 
k Q = —x/2t, and 

. db + 



in the variables x = —x/2, t = t/2. Then 



dt 
dB__ 



1 d 2 B__ Al - 
+ Ab++B 2 _ 



(101) 



2 &£ 2 



The asymptotics of the decaying solutions of this system is well known |L6 
|17| , |18[. The leading terms are 



fi__={) r 1 /2+i, e -ix 2 /2 i - 



(102) 



Here u and £ only depend on k , T, n and B. Substitution of these asymp- 
totic solutions into the differential equations ( |101|) gives 



v = -4u v 

In order to calculate v we use the identity 

OB, 



dx 



2\B__b 



++ 



Let us substitute 



here, 



' — ~ln{l - 2(1 - cos(77)) 7 ^o)(l - 7^o))} 



dx Ant 
From (|0D and (TO) we have 



2LB__6, 



2itt ^o 



2f 



Using fllPl]) , ( |105| ) we get 
1 



2tt 



In {1 - 2(1 - cos(7/)) 7 tf(fco)(l - 7^o))} 



(103) 



(104) 



(105) 



(106) 



(107) 



So we have calculated the asymptotics of b ++ and B , see ( |102| ). We can use 

the differential equations (|101| ) even more. We can improve the asymptotic 
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expression ( p7[ ) for lndet(7 + ^V). Following the steps from pages 455-457 
of the book Q we obtain 

lndet(J + 7 t>) = 

— r dp \x - 2pt\ In (l + j&Op) ( e -™(P- fe o) - l)) + — lnt . (108) 
2% J-oo ^ \ ' ' 2 

I order to get the correlation functions (0), (H) let us calculate 
6 ++ det(J + 7 V") = u r^ 1+i ^ 2 exp(ix 2 /4t) 

• exp y°° dp \x - 2pt\ In (l + 7 0(p) ( e -^gn( P -k ) - l)) | , (109) 

fi__det(J + jV) = vor^ 1 -^ 2 exp(-ix 2 /4t) 

■ exp y°° dp \x - 2pt\ In (l + 7 0(p) (V™^) - l)) } . (110) 

We did not calculate the coefficients uq and w in these expressions. We 
calculated the leading exponential factor and power (of time) corrections. 
Further down, while substituting we shall not pay attention to numerical co- 
efficients. Let us get all notations together, and let us present our asymptotic 
expressions for the correlation functions (|i~lD , (|l~2|), 



G+(x,t) = e ^-*)+i* 2 /4t / dz UoF (z)r^ 1+i ^ 2 e tS ^ ,(111) 

J (z — l) 2 

G^(x,t) = e -m-B)^/« f± VoF ( z ) t -fr-w*)) a e ts(z) (112) 

Here the contour of integration is the unit circle. 

F(z) = 1 + ^ + ^— , (113) 
7 — z 72; — 1 



= -i-ln{l-(2-2-^ 1 ) 7 ^ )(l-7^o))} , (H4) 

tS(z) = — r dp\x- 2pt\ In (l + v&(p) ( z -*zrip-*o) - 1)) , (115) 

2% J-00 ^ ^ ' ' 

7 = 1 + ^ , U = >/a , 7 . W = .(116) 



19 



At the end of section 2 we showed that b ++ cancels the second order pole at 
z = 1. This means that uq will cancel the second order pole at z = 1 in ( |1 1 1| ) . 
The remaining contour integrals in ( |111| ), ( |112j ) can be evaluated by the 
method of steepest descent (see next section). Eqs. ( |111| ), ( |112| ) are integral 
formulae for the large time, long distance asymptotics of correlation functions 
of local fields for the delta-interacting electron gas. These formulae are valid 
for arbitrary finite temperatures. They are the main result of the present 
article. In the next section we shall simplify our asymptotic expressions 
( |1 1 1|) , ( |1 1 2| ) at low temperatures. Let us note that S(z) can be written as 



S(z) = -J_Jp\p-k \ln{ e(fc2+B _, )/T + ~ 1 • (H7) 

7 Low temperature asymptotics 



We obtained the asymptotic formulae (|111|) , ( |112| ). The remaining contour 



integral should be taken by means of the steepest descent (saddle point) 
method. Let us consider here only the leading asymptotic factor. Then we 
can rewrite the simplified expression for the asymptotics as 

G+(x,t) = ^~ B ^e ix2/u n + , (118) 
Gj(x,t) = e-^-^e'^^n. , (119) 

Ot = / dzF(z)t-^ 1±iu{z))2 e ts{z) . (120) 



(121) 



The saddle point equation dS/ dz = can be represented in the form 
dk k poo \z 







1 _|_ J_ e ((fc-fco) 2 +B-^)/T J Q 1 + £ e ((fc+fco) 2 +B-^)/T 



■yz 



k = x/2t, 7 = 1 + e 2B / T . We prove in appendix A that there is only one 
solution of this equation in the interval < z < 1. Considering the pure 
time direction, k = 0, we find the solution 

z = l , (122) 

which leads to 

n± ~ r 1 . (123) 
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This result is valid for all finite temperatures. 

Let us consider the other asymptotic regions, k > 0, x — > oo and t — > oo. 
In these regions the saddle point equation (|121[) can be solved explicitly only 



for small temperatures. There are two solutions z± = ±z c , 

T 3/4 



2^1% 



3/2 



e- fe o/2T (124) 



If we deform the contour of integration from a circle to a line through z + 
parallel to the imaginary axis, we pass the saddle point in the right direction 
(see appendix A). 

In appendix A we obtain the following low temperature expansion for the 
phase tS(z), 

tS(z) = -2k Dt | (l - z 2 c + 1 - z\ . (125) 
Here D is the density @, 

^e"/ T B = , 

D = { v * if (126) 

fZ e (p+\B\)/T B ^ Q _ 



1 



2 

Eq. (125) shows that at low temperatures the relevant parameter for the 



calculation of the asymptotics of Q± is 2koDt = xD rather than t. The pa- 
rameter xD has a simple interpretation. It is the average number of particles 
in the interval [0,a;]. We have to distinguish two limiting cases. 

(a) xD — > 0, the number of electrons in the interval [0,2;] vanishes. In this 
regime the interaction of the electrons is negligible. An electron prop- 
agates freely from to x. Q± cannot be calculated by the method of 
steepest descent. We have to use ( |111|) , (|112| ) to calculate the asymp- 
totics of Gf(x, t) and G7(x, t). Now tS(z) and u(z) tend to zero on the 
contour of integration. Thus Gf(x,t) ~ i~\ e ±^~B)±ix 2 /u ^ w ] 1 j c j lj as 
expected, is the same as for free fermions. The proportionality factor 
is different for G^(x,t) and Gj(x,t). To calculate it from (\1 1 1|) , (|112j) 
we would have to know Uq and Vq. By comparison with free fermions 
(see below) we expect it to be a number for G*(x,t) and proportional 
to e ^- B - k o)/ T for Gf(z,t)- 
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(b) xD ^> 1, the average number of electrons in the interval [0,x] is large. 
Now the interaction becomes important, and we can use the method of 
steepest descent to calculate Q±. This is the most interesting case. In 
the following we will investigate it in detail. 

The contribution of the saddle point to fl±, eq. ( |12U| ), is 

n ± ~ r {1± ' u{z+)) e~ xD ^ , (127) 

where 

D l = |vf e(M+B)/T • ( 129 ) 

Di is the low temperature expression for the density of down-spin electrons, 
eq. (§). 

We further have to take into account the pole of the function F at z = 7 _1 . 
While deforming the contour, we could have crossed it. It turns out that the 
pole contributes to Q±, when the magnetic field B is below a critical positive 
value 

B < B c = kl/A (130) 

(see appendix A). The pole contribution always dominates the contribution 
of the saddle point. It is obtained as 

(3B+M-fc2)/T 

K7- 1 ) = = • (132) 

271 

Finally, we have the following low temperature asymptotic expressions for 
the leading factors of the correlation functions. For 

B<B C : 

G+(x,t) = t -|--(7- 1 ) e «(M-B) e - 2 /4t e -xD i ^ ( 133 ) 
Gf(x,t) = t-|+-(7- 1 ) e -it(^-B) e ~ix' 2 /4t e -xD i ^ (134) 
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where 1/(7 1 ) is given by eq. ( |132| ), and for 
B>B C : 

G+(x,t) = r l--( Z+ ) e i^-B) e - 2 /4t e -^ ^ ( 135 ) 

Gj( X ,t) = t-l+-(^) e -«(M-B) e -- 2 /4t e -x^ ? ( 13 g) 

where v{z + ) is given by eq. (|128|). Recall that we have assumed the aver- 



age number of electrons in the interval [0, x] to be large, xD 3> 1. The 



corresponding expressions for Gf (x, t) follow from equation (|i~3|). 

The leading low temperature expressions for D and _D| agree for positive 
magnetic field. Hence, above the critical field B c the correlation functions 
depend on the magnetic field only through the trivial factors e ±lf ^~ B \ The 
system is saturated. Another way of looking at the condition (|130|) , which 
would be appropriate for experiments with fixed magnetic field B, is the 
following. For positive magnetic field there are two different asymptotic 
regions. One is the space like region, where > AB, the other one is the 
time like region, where k$ < 4B. Correlation functions of local fields have 
asymptotics ( |133| ), ( |134| ) in the space like region and asymptotics ( |135| ), ( |136| ) 
in the time like region. For non-positive magnetic field there is no distinction 
between space like and time like region. 

It is instructive to compare the asymptotic expressions (|133| ) - ( |136| ) 
for the correlation functions with the corresponding asymptotics for free 
fermions. The two-point correlation functions for free fermions have the 
low temperature asymptotics 



tr (e- H °/ T ifo(x,t)i}>\(a,0) 
tr (e- H °/ T ) 

tr (e^°/ T ^{(2,^(0,0) 
tr (e- H °/ T ) 



CVHe^^W 4 ' , (137) 



CJ-V ill ^ B) e- il2/41 , (138) 



as x — > 00, t — > 00. Here Hq is the free Hamiltonian ((|2|) with c = 0). The 
leading low temperature contributions to the constants in ( [L371 ) and ( |138j ) 
are C + = e~ 1 ?/2 v /7r and C_ = e ( - M ~ B ~ ka ^ T e 1 ^ /2y/if. Comparing the free 
fermionic expressions with (|133|) , (|134j) or ( |135| ), (|136|) , respectively, we can 



interpret the factors t ±lv<Kl xDi and t \ ±w< ^ z +">e xC>l as low temperature 
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corrections, which appear due to the interaction. The occurrence of the 
density of down-spin electrons in the exponential factors in eqs. (|133|) - ( |136| ) 
has a natural interpretation. Correlations decay due to interaction. Because 
of the Pauli principle and the locality of the interaction in our specific model, 
eq. (|D, up-spin electrons interact only with down-spin electrons. Therefore 
the correlation length is expected to be a decreasing function of the density 
of down-spin electrons, which diverges as the density of down spin electrons 
goes to zero. We see from ( |133| ) - ( |136| ) that the low temperature expression 
for the correlation length is just and thus meets these expectations. 

In the limit B — > — oo, fi — > — oo, fi — B fixed there are no j-spin electrons 
left in the system, Di — > 0, — > D. This is the free fermion limit. In the 
free fermion limit B < B c , and the asymptotics of G^(x,t) and Gj(x,t) are 
given by eqs. ( |133j ), (|134j) , which then turn into the expressions (|137|) , ( |13S| ) 
for free fermions. 

So we have explicitly evaluated the asymptotics of correlation functions of 
local fields Gf{x,t), (0), flTJ). Our main results are formulae (|1 1 1|) , ( |112| ). 
In this section we considered the saddle point calculation of the remaining 
contour integral, and we discussed the physical meaning of our asymptotics 
in the low temperature regime. 
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Appendix A Steepest descent calculation of 
the leading asymptotic factors 

In this appendix we discuss how to calculate the leading asymptotic contri- 
bution to the integrals in (|111|) , ( |112|) by means of the method of steepest 
descent. The leading contributions determine the exponential decay and the 
power law behavior of the correlation functions (0), fll2|) - We will neglect 
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all sub-leading factors. Then we are left with the calculation of the contour 
integral (OC 



Ot = fdzF(z)t-*Q* b 'We tS W . (A.l) 
The contour of integration is the unit circle. 

F(z) = 1 + ^ + ^— , (A.2) 

7 — 2; 72 — 1 

v(z) = - 7 ^\n{l-(2-z-z~ 1 )^(k )(l-^(k ))\ , (A.3) 

1 /-co f Jk 2 +B-fj,)/T 1 -sign(p-fco) ~) 

SW = -/^^Ip-^lnj e(tI , B .„,) T+7 } ■ (A.4) 

Let us cut the complex plane from —00 to — ^~ x e^ B ^^>l T and from —^e^~ B ^l T 
to zero. The integrand in eq. ( |A.1| ) is analytic in the cut plane except for the 
two simple poles of F(z) at z = j ±l . We may therefore deform the contour 
of integration as long as we never cross the cuts and take into account the 
contributions from the residua, if we cross z = 7 or z = 7 _1 . 

Let us comment on the somewhat unusual fact that t in eq. (|A.1|) appears 
not only in the exponent. Calculating an integral of the form (|A.1|) by the 
method of steepest descent amounts to deforming the contour of integration 
in such a way that we are left with a Laplace type integral. Let us therefore 
consider the large t asymptotics of the integral 

dxt 9(x) e- tf{x) = / dx e -t(m-t9(x)) _ ( A-5 ) 



Here e = \n(t)/t becomes small as t — > 00. We assume that f(x) has a unique 
minimum xq. Then, first calculating the large t asymptotics of the right hand 
side of ( |A.5[ ) by Laplace's method and expanding for small e afterwards, we 
obtain the leading asymptotics 

dxt 9{x) e- tf(x) = ^^t g{xo) -^e- tf{xo) . (A.6) 

This means that the factor f-jC^W) 3 in equation ( |A~ll ) can be treated the 
same way as a t-independent function. Thus, to leading order, a saddle point 
z c contributes to Q± as 

± ^ r |-|(i±iK^)) 2 e -*5(«c) _ (A. 7) 
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The saddle point equation dS/ dz = can be represented in the form 
dk k r°° dkk 



1 + J_ e ((fc-fc ) 2 +B-M)/T J 1 + l e ((k+k )2+B-ti)/T 



(A.8) 
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= x/2t, 7 = 1 + e 2B / T . This is still a transcendental equation, which 
cannot be solved explicitly for z except at low or high temperatures. 

Before turning to the case of low temperatures, let us show that ( A.8|) 
restricted to the positive real axis, z > 0, has a unique solution, which is 
located in the interval < z < 1. Let us define two functions 

r°° dkk /*s 
J±(x)= io 1 + se(***°) a /r • (A ' 9) 
We will further use the abbreviation 

A = - . (A.10) 

7 

Note that A > 0. Using the above conventions the saddle point equation 
( |OD reads 

/M = 5^ = i ■ (A.") 

Now 

f\z) = A{j;(Az)J_(A/z) + (l/z 2 )J + (Az)/(A/z)}/ J 2 _(A/z) , (A.12) 
where the primes denote derivatives with respect to the arguments. Since 

„oo ] ie (k±k ) 2 /T 

J' ± {x) = -j o dk {l + xe{k±ko)2/T)2 < (A.13) 

and J±(x) > for all positive x, we conclude that f(z) is monotonically 
decreasing for positive z. 

We can decompose J-(x) into 

J_(ar) = 2k VTl(x) + J + (x) , (A. 14) 

where 



I(x)=l~dk T ± g? . (A.15) 
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It follows that 

= < x (A16) 

2k VTl(A) + J + (A) " 



since /(a;), J+(a;) > for all positive a;. Thus the saddle point equation ( |A.8| ) 



has no positive real solution z with z > 1. We also learn from ( [A.16Q that 
,2 = 1 solves the saddle point equation for k Q = 0. 

Let us work out the behavior of f(z) for z — > 0+ and 5, /i, T and /c 
fixed. We have the following asymptotic expansions of I(x) and J+{x) for 
large positive x, 



I(x) = >£l + 0(JL) , (A.17) 



J+(x) = - r ' dkke^ k+k ^' T + o(\) . (A.18) 

X Jo \X 2 J 



Hence 



J-(A/z) = (z/A) (k Q Vtff + J™ dk ke^ k+k °y 2 / T ^j + 0(z 2 ) . (A.19) 



We can further show that 



2 



J|m jj + (x)-|(y-Tln(x)-A;o) [=0 , (A.20) 
which implies that for small positive z 

J+(z) = -T\n(z)/2 + O | yj- ln( :) ) . (A.2I) 

Then 

H*) TA / , r ^,^^-1 , „ / */-ln(z) 



2 2 V jo / 

(A.22) 

Eq. ( [Op shows that lim 2 ^ +/(^) = +oo. Since f(l) < 1, eq. flPED , 
we conclude that the saddle point equation (|A.8| ) always has a real positive 
solution in the interval < z < 1. 
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At low temperatures we can solve the saddle point equation ( |A.8| ) ex- 
plicitly. Recall that we are dealing with the gas phase \B\ + \i < 0. In this 
phase 

A- 1 = (i + e -m/T )e {\B\+»)/T (A-23) 

becomes a small parameter as T goes to zero. We will solve the saddle point 
equation ( |A.8|) selfconsistently. Let us assume that 

hr l e~ k ° /T < \z\ < A . (A.24) 

We further assume that k > 0, since the case &o = can be treated for all 
temperatures (cf. eq. (|A.16|) ). Define 

' dkke- {k±ko)2/T . (A.25) 

( |A.24| ) implies that \A/z\ ^> 1, and thus to leading order 

J-(A/z) = ^ ■ (A.26) 
On the other hand \K.z\e k °l T ^> 1, which implies that 

M^) = j- z • (A.27) 

Inserting ( |A.26| ) and ( |A.27| ) into the saddle point equation ( [A. 1 1| ) we obtain 
two solutions z± = ±z c , 



z c = ^I+/I- . (A.28) 
The low temperature asymptotics of I + and /_ are easily calculated. We find 

U - , (A.29) 

J_ = k V^T + I + w fcov^T . (A.30) 

It follows that 

T 3 /4 

z c = - 3/2 e- fc o 2 /^ . (A.31) 
Note that this solution is consistent with ( |A.24j ). 
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Im z 



Re z 



Figure 1: Directions of steepest descent 

In order to determine the directions of steepest descent let us calculate 
S"(z±). First we rewrite S(z), eq. (|A.4|) , as 



1 f» [ 1 + (2:A) -1 e~ (fc+fco) /T 1 + z A~i e -(k-k y/T i 

{z) ~nJo \ 1 + A-i e -( fc + fc o)VT i + A -i e - (fc _ fe0 )VT j • 

(A.32) 

For all z, which satisfy (|A.24f ), we obtain the low temperature expansion 

s < 2 > = Mi- 1 )'^- 1 "-} 

= -2k D | (l - -) z 2 c + 1 - z\ , (A.33) 

where D is the density, eq. (H). Thus 

^) = ±^ , (A.34) 

which means that S"(z+) > 0, S"(z-) < 0. Therefore the path of steepest 
descent through z_ is part of the real axis, whereas the path of steepest 
descent through z + is perpendicular to the real axis. The directions of steep- 
est descent are depicted in figure 1. We conclude from the figure that the 
relevant saddle point for the calculation of our integral ( |A. 1| ) is z + . 

Now the saddle point contribution to the integral ( |A.1|) follows from eq. 
flA~33D . We obtain 



S(z + ) = -2k D(z + - l) 2 w -2k D . (A.35) 
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The exponent v{z), which determines the power law behavior of the correla- 
tion functions, has the following low temperature expansion, 



v \ z ) = 2ti 



± ln(l -(2-z- z-^A-'e-W 7 ) . (A.36) 
Since z c satisfies (|A.24|) and z c < 1, we obtain 



UM = -2^K = - • 



Eqs. ( |A.35|) and (|A.37|) determine the saddle point contribution to fl±, 

Ot ~ r {1±iu{z+)) e~ xD . (A.38) 

Here we suppressed i ,2 (z + ), which is small compared to 1. 

The function F(z), eq. (|A.2|) , has two simple poles at z = 7 ±x , which 
may contribute to the integral (|A.1|) , if they are crossed in the process of 



deformation of the path of integration. Let us investigate this question in 
the low temperature limit. The cases B > 0, B = and B < have to be 
studied separately. 

(a) B > 0: In this case 7 = e 2B//T , 7 -1 = e~ 2B ^ T . Since z + becomes 
small at low temperatures, we may have crossed the pole at 7 -1 . The 
condition for having no contribution from the pole is e~ 2B l T < z + , or, 
using ( |A.31|) , &q < 4B. Let us assume that 



kl > 4B . (A.39) 
Then z = 7 -1 satisfies (|A.24|) , and we can use the low temperature 



expansion ( |A.33| ) to obtain the contribution of the pole, 
Sir 1 ) = -2k D((l- 1 )z 2 c + l- 1 - 1 ) 

w - 2k D(l - -fz 2 c ) w - 2k D . (A.40) 



Here we used (|A.39|) once more to estimate that jz 2 1. We see that 



the contribution of the pole to the exponential decay is exactly the 
same as the contribution of the saddle point. This is no longer true for 
the power law. ( |A.39j ) implies that 7A~ 1 e _fc o/ T <^ 1. Using this fact 



and the low temperature expansion ( |A.36|) we obtain 
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(b) B = 0: In this case 7 = 2, 7 1 = ~ and z + < 7 x . The pole at 7 1 
contributes to the integral (|A.1|) . z = 7 _1 = | satisfies ( |A.24| ). Thus 



5( 7 - 1 ) = - A; 0j D(1 - 2z 2 c ) ^-k D . (A.42) 

Comparing this expression with the saddle point contribution ( |A.35| ) we 
see that the pole now yields the leading contribution to the exponential 
decay. For ^(7 _1 ) we find 

K7- 1 ) = "— = • (A.43) 

Z7T 

(c) B < 0: Now both poles are very close to 1, 7 = 1 + e 2B ^ T , 7 _1 = 
1 — e 2B / T + e 4B//T . The pole at 7 _1 contributes to the integral, and 
z = 7 _1 satisfies ( |A.24|) . Using again (|A.33|) we obtain 



S^ 1 ) = ~ 2k De 2B / T (l - z 2 c ) « -2A; 0J De 2B / T . (A.44) 

Since S < 0, this gives again the leading exponentially decaying con- 
tribution to Q±. ^(7 _1 ) is obtained from (|A.36|) , 

„( 7 -i) = _ _L i n (i + e^/^-ig-fcg/T) w _ £ (A>45) 

Using the explicit formulae for the density D, eq. (ffl), we obtain a uniform 
description of the cases (a), (b) and (c), 



Ot = fdzF(z)r^ l±i ^ 2 e ts W 

„ t -(i±M7- 1 )) exp |_| v /J e (A t +B)/r| ; ( A46 ) 



v{T l ) = ^ • (A.47) 

For a physical interpretation of eq. ( |A.46| ) we look at the formula fl5|) for the 
density of down-spin electrons. Its low temperature expansion is 
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Inserting this expression into ( |A.46| ) we arrive at eq. ( |131| ) of the main text. 
Furthermore, comparing ( |A.48|) and (H) we see that we can identify with 
D, when the magnetic field is positive. Hence we may replace D by D± in eqs. 
( |A.37|) and (|A.38| ), which implies eqs. ( |127| ) and (|128|) in section 7. 



We still have to discuss the range of validity of the above low temperature 
calculations. The low temperature expression for the saddle point z c is suffi- 
ciently precise, as long as z c satisfies (|A.24|) . Yet, there is another condition 
to be satisfied. Note that the phase S(z) in low temperature approxima- 
tion, eq. ( |A.33| ), is multiplied by a factor which vanishes as the temperature 



approaches zero. We have 



tS(z) = -2k Dt (^--j^ + l-.j . (A.49) 

Thus the large parameter, which we need for the saddle point approximation 
to be valid, is xD = 2k^Dt rather than t. 
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